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Abstract 

We investigate analytically and numerically the Multiple Quantum (MQ) NMR dynamics in 
systems of nuclear spins 1/2 coupled by the dipole-dipole interactions in the case of the dipolar 
ordered initial state. We suggest two different methods of MQ NMR. One of them is based on the 
measurement of the dipolar temperature in the quasi-equilibrium state which establishes after the 
time of order T2 after the MQ NMR experiment. The other method uses an additional resonance 
45"0-pulse after the preparation period of the standard MQ NMR experiment in solids. Many- 
spin clusters and correlations are created faster in such experiments than in the usual MQ NMR 
experiments and can be used for the investigation of many-spin dynamics of nuclear spins in solids 
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I. INTRODUCTION 

Multipe-quantum (MQ) NMR spin dynamics in solids [l] is a powerful tool for the in- 
vestigation of structure and dynamical processes in solids, counting the number of spins 
in impurity clusters Q| and the simplification of ordinary NMR spectra [3] . Although MQ 
NMR was successful in a lot of applications and experimental methods of MQ NMR have 
been developed adequately, the theoretical interpretation of many-spin MQ NMR dynam- 
ics is restricted by the phenomenological approach A systematic quantum- mechanical 
approach was developed only for one-dimensional systems in the approximation of 

nearest neighbor dipolar interactions. Up to now the thermodynamic equilibrium density 
matrix in a strong external magnetic field has been considered as the initial condition for 
these experiments and theoretical interpretations. Recently, it has been suggested [3] to 
consider the dipolar ordered state as the initial state for such experiments. It is well known 
that the dipolar ordered state can be prepared using the method of adiabatic demagnetiza- 
tion in a rotating frame (ADRF) 

aa 

or with the Jeener-Broekaert (JB) two-pulse sequence 
s[ [uj] . It has been shown that the spin system contains MQ coherences immediately follow- 
ing the second pulse of the JB pulse sequence . By encoding the coherence numbers in 
an orthogonal basis (the x-basis), it was shown that the dipolar ordered state is a two-spin 
correlated one As a result of using the dipolar ordered initial condition, many-spin 

clusters and correlations appear faster than in the ordinary MQ NMR experiments in solids 
[l| and some peculiarities of MQ dynamics can be investigated with these experiments. Of 
course, it is necessary to make some changes in the scheme of the standard experiment in 
order to obtain non-zero signals of MQ coherences. 

In the present article we consider MQ NMR dynamics when the initial condition is deter- 
mined by the dipolar ordered state. Our motivation of performance of this work is defined 
first of all by that the many-spin correlations are created faster in such experiments and can 
be used for the investigation of many-spin dynamics of nuclear spins in solids. We consider 
peculiarities of MQ NMR in systems prepared in the dipolar ordered states in Section II. 
Two methods of MQ NMR for the systems in the dipolar ordered states are presented in 
Sections III and IV. The method of Section III is based on the measurement of the dipolar 
temperature of the quasi-equilibrium state which establishes after the time t ~ u[ l (uji oc 
is the dipolar local field) after the MQ NMR experiment. The method of Section IV is a 
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simple modification of the usual MQ NMR experiment The additional (7r/4) 2/ -pulse after 
the preparation period allows the observation of a non-zero signal as in the usual MQ NMR 
experiments. Computer simulations of such experiments for linear chains containing up to 
eight spins are presented in Section V. 



II. MQ NMR WITH THE INITIAL DIPOLAR ORDERED STATE 

We consider a system of nuclear spins (s = 1/2) coupled by the dipole-dipole interaction 
(DDI) in a strong external magnetic field. The secular part of the DDI Hamiltonian 8] has 
the following form 



7~Ldz — Djk 

j<k 



(1) 



where = -^(1 — 3cos 2 9jk) is the coupling constant between spins j and k, 7 is the 
gyromagnetic ratio, Tjk is the distance between spins j and k, 9jk is the angle between the 
internuclear vector rjl and the external magnetic field H which is directed along the axis 
z. Ij a {ot = x, y, z) is the projection of the angular spin momentum operator on the axis a; 
1^ and Ij are the raising and lowering operators of spin j. 

The basic scheme of the standard MQ NMR experiments consists of four distinct periods 
of time (Fig. la): preparation (r), evolution (ti), mixing (r) and detection (t 2 ) MQ 
coherences are created by the multipulse sequence consisting of eight-pulse cycles on the 
preparation period jjj . In the rotating reference frame (RRF) the average Hamiltonian 
describing the MQ dynamics at the preparation period can be written as [l| 

n MQ = n^ + n^ 2 \ (2) 

where i?( ±2 ) = —\^2j <k Djklfl k t . Then the evolution period without any pulses follows 
. The transfer of the information about MQ coherences to the observable magnetization 
occurs during the mixing period. The density matrix of the spin system, p(r), at the end of 
the preparation period is 

P (t)=U(t)p(0)U+(t), (3) 

where U(r) = exp(—ir(H^ + H^ 2 ^)) is the evolution operator on the preparation period 
and p(0) is the initial density matrix of the system. On the mixing period the spin system 
is irradiated with sequence of the pulses shifted on a ir/2— phase regarding the pulses of the 
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preparation period. As a result the average Hamiltonian describing evolution of the spin 
system on the mixing period changes a sign on opposite to the sign of the Hamiltonian (2), 
and the evolution operator , U(t), is replaced by the operator U + (r). Usually the thermo- 
dynamical equilibrium density matrix is used as the initial one for MQ NMR experiments. 
Here we consider MQ NMR dynamics with the initial dipolar ordered state when the Hamil- 
tonian of the system is determined by Eq. (1). Then the equilibrium state can be described 
as 

p(0) = 1 ex P (-(3H dz ) « (1 - (3H dz ) /2 N , (4) 

where j3 is the inverse spin temperature, the partition function Z = Tr{exp(—j3Hd z )} and 
N is the number of spins in the system. The high temperature approximation is taken into 

3e reached by using 



{ |9| or by appl ying 



account in Eq. (4). The dipolar ordered state of the spin system can 
the method of adiabatic demagnetization in a rotating frame (ADRF) 
a pair of phase-shifted radiofrequency (rf)-pulses (the Jeener-Brokaert (JB) method) [a. Il0]. 
It is evident that the unit operator in Eq. (4) is not significant for the time evolution of the 
density matrix. For simplicity we will take 

p(o) = n dz (5) 

as the initial condition. The use of the dipolar ordered state as the initial condition in MQ 
NMR experiments leads to the emergence of spin clusters with a greater number of spins 
faster than in the standard method. It is well known |8j that NMR methods in solids use 
both the Zeeman equilibrium state and the dipolar ordered one as the initial conditions. 
Here the developed methods together with the standard one Q yield analogous possibilities 
for MQ NMR spectroscopy. Below we describe two different methods of MQ NMR of the 
dipolar ordered systems. 

III. THE DIPOLAR TEMPERATURE AS A SOURCE OF INFORMATION 
ABOUT MQ NMR 

Starting with the initial condition (5), the density matrix, p(t), after the three periods of 
the standard MQ NMR experiment (Fig. la) can be written as 

p(t) = U + (r)e- iStlh p(T)e l5tlh U(r), (6) 
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where p(r) is the density matrix at the end of the preparation period according to Eq. (3) 
and t — 2t + ti, 6 is the frequency offset on the evolution period of the duration U which 
is a result of applying the time proportional phase incrementation (TPPI) method The 
last unitary transformation in (6) with operator U(r) describes the mixing period. Then 
the evolution of the system in the RRF is governed by the Hamiltonian H dz of Eq. (1) and 
after the time t > lo^ oc (u>f oc = Tr{7i 2 dz ] /Tr{I 2 }) the system achieves the two-temperature 
quasi-equilibrium state with the density matrix p eq , 

p eq = au I z + pH dz , (7) 

where the inverse temperatures, a and (3, can be found from the conservation laws: 

_ Tr{p(t)I z } _ Tr{p{t)H dz } 

u,Tr{P z } ' P Tr{Hl} ' [ *> 

Applying the unitary transformation V = e ~' l ^ Iz e~' lnIx to the expression p(t)I z one can 
obtain from Eq. (8) that a = 0. It means that the quasi-equilibrium state of Eq. (7) is the 
dipolar ordered state. It is convenient to expand the density matrix, p(r), at the end of the 
preparation period of the MQ NMR experiment (Fig. la) with the initial condition (5) as 

1 

p(r) = J2Pn(r), (9) 

n 

where the term p n (r) is responsible for the MQ coherence of the n-th order. One can find 

e- m ^p n (r)e m ^ = e- m5t p n (r). (10) 
By using Eqs. (8)-(10) one can rewrite the temperature (3 as 

p = J2e~ m5t J n (r), (11) 

n 

where the intensities of MQ coherences 



Tr {p n (r)p- n ( r)} 



J n {r) = ±,W ^ ( Z^ — ■ (12) 



The temperature (3 can be measured using the magnetization component which is in 
phase with the exciting rf-pulse Q|. Different frequency components of the temperature (3 
yield MQ coherences of different orders. 
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IV. THE MODIFICATION OF THE STANDARD MQ NMR EXPERIMENT IN 
SOLIDS 



In the previous section we showed that at the initial condition of Eq. (5) the standard 
scheme of MQ NMR experiments does not lead to any observable magnetization. In order 
to overcome this problem we introduce a 0^-pulse turning spins around the axis y by the 
angle <fi after the preparation period (Fig. lb). Then starting with p(0) = Hd z , after the 
mixing period (Fig. lb) the longitudinal magnetization, M z (t) , is 

M z (t) = Tr{p(t)I z } (13) 

where the density matrix after the three periods of the modified MQ NMR pulse sequence 
(Fig. lb) can be written as 

p(t) = U + (T)e- i5IA p^(r)e i5hh U(T) (14) 
and p^r) is the density matrix just after the y -pulse (Fig. lb) 



p^r) = e-^f/(r)p(0)f/ + (r)e^. (15) 

It is worth to notice that here the TPPI method [l| is applied in the mixing period in 
contrast to the method of Section III when it was applied in the preparation period of the 
MQ NMR experiment. Using Eqs. (14) and (15) and the initial condition (5) it is convenient 
to present the formula (13) for the longitudinal magnetization, M z (t) , as follows 

M z (t) = Tr {e-^e-^U^H^U+iry^e^pMQiT)} (16) 

where 

Pmq (t) = U(t)I z U + (t), (17) 

coincides with the density matrix at the end of the preparation period of the standard 
MQ NMR experiment with the thermodynamical equilibrium density matrix as the initial 
condition [1]. The density matrix p MQ (r) can be represented in the following form [l] 

P MQ {r) = Y.Pn Q ^) (18) 

n 
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where again the term Pn®{j) is responsible for the MQ coherence of the n-th order and 
satisfies to the relationship of Eq. (10). By using Eqs. (10), (17), and (18) one can rewrite 
the expression for the observable signal in terms of the intensities of MQ coherences 

M z (t) = -i £ e- in5t J n (r) = i £ e m5t r n (r) = i £ e"*** J* n (r), (19) 

n n n 

where 

J n (r)=tTr{p,(r)p^(T)} (20) 

and the longitudinal magnetization, M z (t), is always real. Thus, one can find from Eq.(19) 
that 

J n (r) = -r n (r). (21) 

One can notice that 

I Z (S) = -U-5) (22) 

by applying a 7r-rotation around the y-axis to all operators in Eq. (13). It follows immedi- 
ately from Eqs. (19) and (22) that 

Ju{t) + J- n (r) = 0. (23) 

It is evident that with the method of Section III J n (r) = J_ n (r), and the longitudinal 
magnetization, M z (t), is determined as 

M z (t) = ^2cos(n5t)J n (T). (24) 

n 

By using the modified pulse sequence (Fig. lb) the longitudinal magnetization can be 
expressed as 

M z (t) = ~Y^ sm(n5t)J n (T). (25) 

n 

Thus the phases of the MQ NMR coherences in the modified MQ NMR pulse excitation 
scheme are shifted by re/ 2 from the phases in standard MQ NMR [l]. It is clear from Eq. 
(23) that the intensity of the MQ coherence of the zeroth order, Jo(t), is equal to zero in 
the considered experiments. 

It is well-known that in the usual MQ NMR experiments the sum of the intensities of all 
MQ coherences does not depend on time [l^] . This is also right for the method described in 
Section III. Here this law has a specific form. According to Eq. (23) the sum of intensities 
of orders n and — n is equal to zero for all n. Thus, the sum of the intensities of all MQ 
coherences equals zero. 



V. THE NUMERICAL ANALYSIS OF THE TIME EVOLUTION OF MQ COHER- 
ENCES 



We restrict ourself to numerical simulations of MQ NMR dynamics of one-dimensional 
systems. For example, quas i-one-dimensional hydroxyl proton chains in calcium hydroxya- 
patite C a§{OH)(PO/±)z [lj] and fluorine chains in calcium fluorapatiteCa5F(P04)3 are 
suitable objects to study MQ dynamics in dipolar ordered systems. The numerical calcula- 
tions were performed with the methods of Sections III, IV for MQ NMR dynamics of linear 
chains of 6 and 8 spins. The DDI coupling constant of the nearest neighbors is chosen to be 
D = Is -1 . Then the coupling constant of spins j and k is D/\j — A;| 3 . In order to compare 
the results of the numerical simulations with the analogous ones for the ordinary MQ NMR 
dynamics we introduce the normalized intensities of MQ coherences for the method of Sec- 
tion IV, J n (r)/{Tr(/^)Tr(ifJ 2 )} 1 / 2 . The dependence of the intensities of MQ coherences on 
the dimensionless time, t = Dt , of the preparation period in spin chains containing six and 
eight spins is presented in Figs. 2 and 3 for the both methods of Sections III and Sections 
IV. One can compare the intensities of MQ coherences in the suggested experiments with 
the standard ones (see Figs. 2 and 3). It is evident that the suggested methods can be 
considered as a useful addition to the standard MQ NMR methods. Comparison of the 
suggested methods of MQ NMR of the systems in the dipolar ordered state is given in Fig. 
2 for MQ coherences of the second and fourth orders. One can see (Figs. 2a, 2b) that the 
method of Section III is more preferable for experimental realizations in some cases. 

Fig. 2a demonstrates that the method of Section III yields the intensities of MQ coher- 
ences of the fourth order which several times higher than the analogous coherence in the 
standard MQ NMR experiment. It is clear from the inset of Fig. 3a that MQ coherence of 
the sixth order, obtained with the method of Section IV, appears little earlier than in the 
usual MQ NMR in a linear chain of six spins. At the same time, the intensity of this MQ 
coherence obtained with the method of Section III is equal zero. In fact, the both methods 
supplement each other at the study of MQ NMR dynamics. A tendency of the faster growth 
of MQ coherences of high orders takes place for the linear chain containing eight spins (Fig. 
3b). This peculiarity is connected with the initial dipolar ordered state, p(0) = H dz . As a 
result, many-spin clusters and connected with them MQ coherences appear faster than in 
the standard MQ NMR with the initial condition, p(0) = I z . The numerical calculations 
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confirm also the results obtained in the previous section. In particular, the growth of MQ 
coherences (for the method of Section IV) occurs in accordance with Eq. (23) and Jo(t) = 0. 

One can see from Figs. 2 and 3 and Eq. (23) that the observable intensities can be 
negative in contrast to the ordinary MQ NMR experiments at high temperatures 
At the same time, it was shown [15j that the intensities of MQ coherences can be negative in 
the standard NMR experiments at low temperatures. It is a terminological problem only. In 
fact, in MQ NMR experiments the observable is the longitudinal magnetization modulated 
by rf pulses. The distinct frequency components of the magnetization can have an arbitrary 
sign (la ]. 

VI. CONCLUSIONS 

The MQ NMR methods for the detection of MQ coherences starting from the dipolar 
ordered state are proposed. The first method is based on the quasi-equilibrium state which 
establishes after the time t > oj^ oc after the MQ NMR experiment. This state is the dipolar 
ordered one and its temperature is a source of information about MQ NMR dynamics. In 
order to observe a non-zero signal of MQ coherences in the second method an additional 
resonance (7r/4) y -pulse should be applied after the preparation period of the standard MQ 
NMR experiment in solids. Investigations of MQ NMR dynamics in the dipolar ordered 
states can be considered as a supplementary method which complements the usual NMR 
in order to study structures and dynamical processes in solids. Many-spin clusters and 
many-spin correlations are created faster in such experiments than in the usual MQ NMR 
with the initial condition without any correlation between the spins. This paper focuses 
on simple one-dimensional examples but the physical picture obtained here is not limited 
to the performed simulations and experiments and opens new possibilities for the study of 
many-spin systems. 
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Caption figures 

Fig. 1 The basic schemes for MQ NMR: a) the scheme of the standard MQ NMR 
experiment; b) the scheme of the modified MQ NMR experiment. 

Fig. 2 The time dependence of the intensities of the MQ coherences of the second order 
in a linear chain of six spins coupled by the DDI for p(0) = H dz ; the intensities J_ 2 (solid) 
and J + 2 (dash) obtained by the method of Section IV at = 7r/4. The intensity J2 (dot) 
is obtained by the method of Section III. (b) The time dependence of the intensities of the 
fourth order MQ coherences in a linear chain of six spins coupled by the DDI : the intensity 
of the MQ coherence of the fourth order (solid) for p(0) = I z ; J_ 4 (dash) and J +4 (dot) 
for p(0) = H dz obtained by the method of Section IV at = 7r/4, and the intensity J 4 
(dash-dot) is obtained by the method of Section III. 
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Fig. 3 The time dependence of the intensities of the sixth order MQ coherences in a linear 
chain of six spins coupled by the DDI: the intensity of MQ coherence of the sixth order (solid) 
for p (0) = I z ; J_ 6 (dash) for p (0) = H dz obtained by the method of Section IV with using 
(7r/4) J/ -pulse and Jg(dot) obtained by the method of Section III. (b) The time dependence of 
the intensities of the sixth order MQ coherences in a linear chain of eight spins coupled by 
the DDI : the intensity of MQ coherence of the sixth order (solid) for p (0) = I z ; J_6 (dash) 
for p(0) = H dz obtained by the method of Section IV with (7r/4) y -pulse; J 6 (dot-dash) and 
J 8 (dot) obtained by the method of Section III. The insets show that the MQ coherences in 
the dipolar ordered state obtained by the method of Section IV with applied (7r/4)j / -pulse 
appears little earlier than in the usual MQ NMR and in the method of Section III. 
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